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Abstract 

The higher twist contributions to the deeply inelastic structure functions ^^(xjQ^) and 
(s^) Q^) foi^ larger values of the Bjorken variable x are extracted extrapolating the twist- 
2 contributions measured in the large W region to the region 4 GeV^ <W^ < 12.5 GeV^ 
applying target mass corrections. We compare the results for the NLO, NNLO and N^LO 
analyzes and include also the large x at N^LO to the Wilson coefficients. A gradual lowering 
of the higher twist contributions going from NLO to N^LO is observed, which stresses the 
importance of higher order corrections. 



Deeply inelastic structure functions contain higher twist corrections [1] both in the region of 
large and small values of the Bjorken variable x [2-5]. While the leading twist sector both 
for unpolarized and polarized deeply inelastic scattering is well explored within perturbative 
Quantum Chromodynamics (QCD) up to the level of 3-loop, resp. 2-loop, corrections [6,7], very 
little is known on the scaling violations of dynamical next-to-leading twist correlation functions 
and the associated Wilson coefficients [1] , even on the leading order level. ^ In many experimental 
and phenomenological analyzes, cf. [2,3], higher twist contributions are parameterized by an 
'Ansatz' [2], which is fitted accordingly. Within QCD this ad-hoc treatment usually cannot 
be justified, performing at the same time a higher order analysis for the leading twist terms. 
Since neither the corresponding higher twist anomalous dimensions nor Wilson coefficients were 
calculated, the data analysis has to be limited in the first place to the kinematic domain in which 
higher twist terms can be safely disregarded. 

In the case of fiavor non-singlet combinations of structure functions this is widely the case 
in the region 

g2 > 4GeV^,W^ > 12.5 GeV^ as detailed analyzes of the large region show, 
cf. [12]. In this region one may perform a three-loop QCD analysis, which requires the 0{al) 
Wilson coefficients [14] and the 3-loop anomalous dimensions [6]. The analysis can even be 
extended effectively to 4-loop order, since the dominant contribution there is implied by the 
3-loop Wilson coefficient [7], parameterizing the yet unknown 4-loop anomalous dimension with 
a ± 100 % error added to an estimate of this quantity formed as Fade- approximation out of 
the lower order terms. A comparison with the 2nd moment of the 4-loop anomalous dimension 
calculated in [15] showed [12] that the agreement is better than 20 %, which underlines that the 
above approximation may be possible. We limit the QCD-analysis of the twist-2 contributions 
to this representation since neither as{iJ>^) nor the splitting and coefficient functions are known 
beyond this level. 

The evolution equations are solved in Mellin-iV space. The non-singlet structure function at 
the starting scale of the evolution, Qq, is given by 

oo 
fe=0 

with C^^{N) the expansion terms of the non-singlet Wilson coefficient with Co{N) — 1, 
O'siQ'^) = «s(Q^)/(47r) and /^''^''^'^(iV, Q^) the corresponding combination of quark distributions, 
cf. [12]. Here we identify both the renormalization and factorization scale with Q^. Beyond 
0{al) dominant large x contributions to the Wilson coefficient were calculated in [13]. 
The evolution equation for the quark densities to 4-loop order reads : 



\aoJ 




""^In the small x region saturation eorreetions, see [4,5], are considered which also belong to the class of 
higher twist corrections, applying different formalisms. An important question to be addressed here concerns 
sub-leading corrections, which may be large, cf. [8]. This has been confirmed in a fixed order calculation to three 
loops in [9] . Recent investigations of the small-a; behaviour of deeply inelastic structure functions account for 
these efi'ects [10,11]. 
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Here Pj^{N) denotes the {k + l)-loop anomalous dimension and Pk are the expansion coefficients 
of the QCD /3-function, of. [16], with 



(3) 



k=0 



Eq. (3) is solved perturbatively to 4-loop order in the MS scheme [17] observing the flavor 
matching conditions for the renormalization scale /i at the thresholds rric = 1.5 GeV and mb = 
4.5 GeV, respectively, to be able to compare to other measurements of as{M'^), resp. Aq^^j-,. To 
perform the data analysis the expression for the structure functions, (1), is transformed back to 
x-space by a numeric contour integral around the singularities of the problem in the complex 
A'"-plane. 

In the analysis mass corrections have to be accounted for. These are the target mass [18] 
and heavy flavor corrections [19-21]. While the former are signiflcant, the latter ones amount 
only 1-2 % at NLO and are expected to be even smaller in the yet unknown higher orders^ 
in the flavor non-singlet case. In the evolution equations (1,2) the anomalous dimensions and 
coefficient functions are represented in A^-space [6,7,20,24,25]. Here we applied simplifications 
due to algebraic [26] and structural relations [27] between harmonic sums. Under the conditions 
mentioned above we perform the twist-2 analysis from leading order (LO) to 4-loop order (N^LO) 
fitting the non-singlet ^^''^(a;, Q^) world data, cf. [12], using MINUIT [28]. We then extrapolate the 
results to the region 4GeV^ < W'^ < 12.5 GeV^ and determine effective higher twist coefficients 
Cut{x, Q^) given by 



+ 



g2[i GeV\ 



(4) 



Here Otmc[ ] denotes the operator of target mass corrections. 

QCD corrections beyond N^LO are known in form of the dominant large-x contributions to 
the QCD-Wilson coefficients [13]. Since these corrections do quantitatively only apply in the 
range of large x we do not use them in the twist-2 QGD-fit, because here the data are mainly 
situated at lower values of x where beyond 4-loop order other contributions to the Wilson 
coefficients, which are not calculated yet, are as important. Furthermore, the 4-loop anomalous 
dimensions were not calculated yet. The leading large x contributions are given in terms of 
harmonic sums of the type 5*1^1^, ..^i(A^) which obey a determinant representation [29] in single 



^First contributions which arc relevant for the O(a^) heavy flavor contributions to F2{x, Q^) were calculated 
in [22] for the region m? . Under this kinematic condition the corresponding corrections to F]^{x, Q^) were 

calculated in [23]. 
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harmonic sums Si{N), i.e. they are polynomials of single harmonic sums.^ One may calculate 
these sums recursively, [29] . The n-fold sum reduces to the {n — A;)-fold sums by 



(5) 



To obtain the large x behaviour we retain the terms 



Si{N) oc ln(A^)+7£ 
Si{N) (X 0, ^>2, 



(6) 
(7) 



as |A^| — > oo, with 7^ the Euler-Mascheroni number, and Q the Riemann ^-function at integer 
values. We agree with the numerical parameterization given in Table 1, [10a]. We take into 
account the N^LO terms in this approximation, which are added to the twist-2 fit results in 
N^LO extrapolating to the region of lower values of W'^. 

The effective higher twist distribution functions C^!^{x) extracted arc shown in Figures 1 and 
2 from NLO to N^'LO. Here we averaged over the values in within the x-bins. ^ The leading 
twist terms are those given in [12], with the values of A^^^d = 265 ± 27, 226 ± 25, 234 ± 26 MeV, 
resp. in NLO, NNLO, and N^LO. Both for the proton and deuteron data Cht(2^) grows towards 
large values of x, and takes values ~ 1 around x = 0.6. The inclusion of higher order corrections 
reduces Cut{x) to lower values with a gradually smaller difference order by order. Yet for the 
highest bins, x > 0.8, the effect of the large x resummation terms is important. Earlier higher 
twist analyzes [3] limited to the next-to-leading order corrections are thus corrected by factors of 
2 and larger at large x to lower values. In the present analysis we limited the investigation to the 
inclusion of the large x terms in N'^LO which are still in the vicinity of a nearly complete QCD 
analysis as outlined above. The present description is likely to be final for values of a; < 0.8. 
Beyond this range there are only few data. More data in this interesting region would be 
welcome and can be obtained at planned high-luminosity coUiders such as ETC [31]. Data from 
measurements at JLAB cannot be included into the present analysis, since the kinematic region 
covered currently averages over the resonances applying duality, which may lead to different and 
probably lower higher twist contributions, cf. [32]. 

In the present analysis we extracted the large x dynamical higher twist contributions to the 
structure functions in a model-independent way. It would be interesting to compare moments 
of the term CHT(a^) to lattice results, which allow to simulate the moments of the corresponding 
higher twist correlation functions, in the future. 
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^Related representations were given in [30]. 

^ We took the opportunity to re-bin the data at very large x, if compared to [12] , to optimize w.r.t. experimental 
errors. 
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Figure 1: Comparison of the higher twist coefficient Cht(^) in the large x region for the proton 
data as function of x in a NLO (dotted hue), NNLO (dashes hue), N^LO analysis (dash-dotted 
line) and adding the large x terms in 0[a^^ for the non-singlet QCD Wilson coefficient (full 
line). Some bin centers are slightly shifted for better visibility. 
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Figure 2: The coefficient Cht(3;) for the deuteron data. The curves have the same meaning as 
in Figure 1. 
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